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Abstract 

Sr2Ru04 is at present the best candidate for being a superconducting analogue of 
the triplet superfluidity in 3 He. This material is a good (albeit correlated) Fermi liq- 
uid in the normal state and an exotic superconductor below T c . The mechanism of 
superconductivity and symmetry of the order parameter are the main puzzling issues of 
on-going research. Here we present the results of our search for a viable description of the 
superconducting state realised in this material. Our calculations are based on a three- 
dimensional effective three-band model with a realistic band structure. We have found a 
state with non-zero order parameter on each of the three sheets of the Fermi surface. The 
corresponding gap in the quasi-particle spectrum has line or point nodes on the a and 
/3 sheets and is complex with no nodes on the 7 sheet. This state describes remarkably 
well a number of existing experiments including power low temperature dependence of 
the specific heat, penetration depth, thermal conductivity etc. The stability of the state 
with respect to disorder and different interaction parameters are also analyzed briefly. 

Introduction. The discovery of superconductivity in S^RuCU |ll EJ has generated renewed 
interest in this material ||. Originally the main motivation in Ref. Ill] was the close similarity 
of its crystal structure with that of high temperature superconducting oxide L^-xCa^CuCU. 
However, unlike cuprates, strontium ruthenate is metallic without doping and is not readily 
superconducting. Indeed, very pure single crystals have superconducting transition temper- 
ature T c = 1.5K, rather low on the scale of the high T c cuprates. 

In the normal state the electrons in strontium ruthenate form a well behaved, albeit aniso- 
tropic and correlated, Fermi liquid. In particular the resistivity shows a typical quadratic 
dependence on temperature pi(T) — pi + A{T 2 at low temperatures. Here pi is the 'residual' 
resistivity measured at temperature just above T c and i denotes the direction of the current 
(i = ab for in plane and i = c for c-axis resistivity). Electronic transport is very anisotropic 
with the ratio of p c /p a b exceeding 500. Normal state electronic specific heat shows a typical 
linear dependence on temperature, C e = jT, with a high value of 7 = 37.5 mJ/K 2 mol. The 
appreciable enhancement (factor 3-5) of 7 over the band structure calculations indicates the 
existence of strong carrier correlations in the system. The same conclusion can also be drawn 
from the enhanced value of the ratio between Pauli spin-susceptibility and 7 (Wilson ratio) 
and the so called Kadowaki- Woods ratio of the coefficient A a i> to 7. These values exceeds 
those for non-correlated systems. 

Consistent with the above picture is the fact that first-principles calculations based on the 
Local Density Approximation (LDA) give a good qualitative account of the electronic struc- 
ture 0]. The energy spectrum around the Fermi level of strontium ruthenate is dominated 
by electrons occupying t2 g orbitals {d xz ,d yz ,d xy ) of Ru i+ . The Fermi surface consists of three 
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cylindrical sheets open in the z-direction. Those called a and (3 stem from hybridised d xz and 
d yz orbitals, while the 7 sheet is mainly of d xy character. The a sheet is a hole Fermi surface. 
It is centered at the point X of the bet Brillouin zone and possesses two fold symmetry. The 
calculated shape of the Fermi surfaces agrees nicely with that measured via de Haas-van 
Alphen effect Q. 

While the normal state, as shown above, is quite typical of a good metal, the supercon- 
ducting state is very exotic. The superconducting transition temperature T c is strongly sup- 
pressed by even small amount of impurities ||. The NMR measured relaxation rate shows 
no Hebel-Slichter peak, and the /iSR experiments indicate the appearance of spontaneous 
magnetic fields at T < T c . Together with a temperature independent Knight shift, all that 
points towards very unconventional, presumably spin-triplet, and superconductivity with a 
time reversal symmetry breaking pairing state. 

Because the electrons forming Cooper pairs in superconductors are fermions, the pair 
wave function has to be odd with respect to interchange. It consists of spin and orbital 
components. If the total spin of a pair is zero then such a pairing state is called spin-singlet. 
The corresponding spin-wave function is odd so the orbital part has to be even. In Galilean 
invariant systems they correspond to even values of the orbital quantum number I — 0,2 . . . 
and the corresponding superconducting states are referred to as having s—,d—, . . . wave 
symmetry. The s— wave symmetry is realised in conventional superconductors, while d— wave 
is the case in high temperature superconductors. The electrons in strontium ruthenate most 
probably pair in relative spin triplet state. Their spin wave function is even and thus the 
orbital one has to be odd. Even though a crystal has only discreet rotational symmetry, the 
corresponding symmetries of the superconducting state are still called p— , f— , . . . wave. 




Figure 1: Calculated tight binding Fermi surface of Sr2Ru04. 

As well known, a relatively weak dependence of T c on the concentration of impurities and 
an exponential dependence of the electron specific heat on temperature point to spin singlet 
s-wave superconductivity. On the other hand a power law temperature dependence of specific 
heat is characteristic of a superconducting gap which vanishes at points or along lines on the 
Fermi surface, and a dramatic decrease of T c with impurity concentration is a signal of I > 
'wave' paring. 
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Table 1: Irreducible representations in a tetragonal crystal. The symbols X, Y Z represent 
any functions transforming as x, y and z under crystal point group operations, while I 
represents any function which is invariant under all point group symmetries. 
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No 
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E u 
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The puzzling behaviour of S^RuCU in the superconducting state is connected with the 
facts that NMR, /j,SK and related experiments indicate the realisation of the A(k) = e z (k x ± 
iky) state which does not vanish at the Fermi level and that the power law temperature 
dependence of the specific heat |Q, penetration depth || or thermal conductivity |J which 
require the gap to vanish along lines of the Fermi surface. The point is that out of all the 
symmetry distinct states of a bet crystal [^) none of the odd-parity states have to be time 
reversal symmetry breaking and also possess gap nodes at the same time (c.f. table I). The 
various states, formally fulfilling both requirements proposed in literature [l2| can be 
shown to be the sum of symmetry allowed states with distinct transition temperatures. The 
single superconducting transition observed in all studied samples rules out all these proposals 
as the single transition could only be a result of accidental degeneracy. Instead of relying on a 
symmetry arguments alone, we present a methodology which is based on explicit construction 
of an effective pairing interaction. 

Below, we shall present in Sec. 2 the model and our approach. In Sec. 3 we present the 
results of calculations, and finally we end with our conclusions. 



The model and the approach. Our model of superconductivity in strontium ruthenate 
is motivated by the experimental facts, summarised above. The electronic structure features 
three Fermi surfaces, and all of them are gaped to ensure the vanishing of the specific heat at 
T = OK. Further, there exists a single superconducting transition and the gap has to both 
break the time reversal symmetry and vanish on the lines of the Fermi surface. We shall 
fulfill all these constraints in the context of a realistic three dimensional band structure with 
parameters fitted to the known details of the Fermi surface, and assumed effective attractive 
interactions between electrons occupying various orbitals. 

We thus take the following simple multi-band attractive Hubbard Hamiltonian: 

H = ^ ^ ((Sfn H)$ij$mm' ^mm' (ij)) ^ima^jm' a 

ijmm' ,cr 

~ 2 ^ ] U mrn i {ij)nimirnjm' a 1 (1) 
ijmm' a a' 
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Figure 2: The relative positions of Ru ions in the tetragonal Sr2RuOn lattice. The heavy 
lines symbolise the main interactions. Interplane Uu for the electrons occupying d xy orbitals 
at the neighbouring in planes ions and Uj_ for out of plane neighbours and d xz and d yz 
orbitals. 

Table 2: The values of parameters used in calculations. Here t, t' are in-plane hopping 
integrals between c (d xy ) orbitals, t ax ,t ox are in-plane hopping integrals between a (d xz ) 
orbitals along x and y directions, t a i, is the in-plane hopping between a and b orbitals 
along x + y, and t^y^ and tj_ are out of plane hopping integrals along the bet body-centre 
vector from a(b) to c and from a(b) to b(a), respectively. 
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as our model which is to describe superconductivity in Sr2Ru04. Here m and mf refer to the 
three Ruthenium t2 g orbitals to be denoted a = d xz , b = d yz and c = d xy in the following, i 
and j label the sites of a body centered tetragonal lattice. The hopping integrals t mm ' (ij) and 
site energies e m were fitted to reproduce the experimentally determined three-dimensional 
Fermi Surface |l^]. We found that the set t mm i shown in Table 2 gave a good account of 
the data. 

In choosing the interaction parameters U^^i we adopted a frankly semi phenomenological 
approach. Namely, we eschewed any effort to derive these from an assumed physical mecha- 
nism of pairing and endeavoured to find a minimal, hopefully unique, set of parameters which 
describes the available data. The point of such strategy is that at the end we can claim to 
have identified the position and orbital dependence of the interaction and therefore provided 
guidance for constructing microscopic models for specific mechanisms. 

In the above spirit we consider two sets of interaction constants: [/" / for nearest neigh- 
bours in the plane and (7 7 ^ m / for nearest-neighbour Ru-atoms in the adjacent plane as indi- 
cated in Fig. 2. To limit further the parameter space we need to explore, we assume that 
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Figure 3: Lowest energy quasiparticle eigenvalues, £^(k) on the 7 (a) and /3 (b) Fermi 
surface sheets; (a) polar plot of the 7 sheet in the plane k z = 0, E 7 , m0 i 0<-f) = 0.22meV, 
E lt mini}^F) = 0.056meV (b) vertical cross-section of the cylindrical /3 sheet in the plane 
k x = Ep tTnax (\ip) = 0.32meV. A non-zero /—wave order parameter (dotted line) lifts 
the p-wave line nodes (solid lines). 



Uf = U\\ , UT = Uf = U± and contemplate only 

/ u u u \ / U± U± u' \ 

Ul m , = 11 u u and Ut m . = \ U x U± u' \ (2) 

y u u [/" / y u u u J 

In fact, we take u = u' = for most of our calculations and we use only one in plane, 
U\\, and one out of plane Ux, parameters to fit the experimental data on the specific heat, 
Cv(T), superfiuid density n s (T) and the thermal conductivity kt(T). The purpose of the few 
calculations with u 7^ 0, u 7^ was only to investigate the stability of the U«, Ux results to 
variations in U^, and U^ m , . 

Within the above model, we solved the Bogolubov-de Gennes equations: 

V- ( E" ~ H m>m f (ij) A™' m ,(ij) \f v? jma , \ 
^ \ A w(y) E" + H mm ,{ij) ) \ v» m , a , J 

together with the self-consistency condition 

A mm ' ^mm'(y)Xmm'(y)i Xmm' ) ^ ^ u ima v jm' a' )), (4) 

which follow from Eq. (1) on making the usual BCS-like mean field approximation Q]. Since 
H mm i(ij) and A^ m ,(ij) depend only on the difference Ri — Rj the solution of Eq. (3) is 
rendered tractable by taking its lattice Fourier transforms and, thus transforming the problem 
into that of a 12 x 12 matrix, eigenvalue problem at each k in the appropriate Brillouin Zone. 
For u = u' = the general structure of Al: , (k) turns out to be of the form 

A cc (fc) = A* c sin k x + A^ c sin k y 

A™™/ (k) = AfL™/ sin cos — cos — + A^ , sin — sin — sin — — 

iiliil \ / mm c\ c\ c\ 1 mm OOO 

( a X • ky . y ky ky \ k Z C , , 

+ A™.™/ sin — cos — + A y , cos — sin — cos (5) 

\ 2 2 mm 2 2/2 w 
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Figure 4: Calculated specific heat, 17, as a function of temperature, T, compared to the ex- 
perimental data of NishiZaki et al. M . The dotted-line includes the /-wave order parameter, 
the solid line is without /-wave. 



for m,m' = a or b. 

Our self-consistency procedure always converged on the amplitudes A^ c , A^ c , Aj^ m /(T), 

A im'( T )' A m,m'( T ) and A ™ m ' ( T ) to an accuracy higher than 10" 4 %. The principle results 
of such calculations are these amplitudes and the corresponding quasiparticle energy eigen- 
values E v {k). These eigenvalues are shown in Fig. 3, for the interaction values Uu, U± given 
in Table 2, and with the constraint that A^j.nfT) = 0. Evidently from the figure, on the 7 
sheet there is an absolute gap below -E 7jm i n (fc_F), in the quasi-particle spectrum, while the 
f3 sheet is gap less with a line of nodes in the gap functions at k z = ±— This dramatically 
different gap structure and symmetry on different sheets of the Fermi Surface is the striking 
new results of an interaction matrix U™ n i (ij) which couples electrons in different Ru-planes. 

From the results in Fig. 3, one is encouraged to investigate this form of interaction further 
because, on the one hand, the line of zero gap on the (3 sheet explains the power-law behaviour 
of various thermodynamic quantities and, on the other, on the fully gapped 7 sheet A 77 (fc) = 
A cc (fc) = AJ c (sin k x a + i sin k y b) which implies the broken time reversal symmetry demanded 
by a number of other experiments. Indeed, for the above pair of interactions, U± and Uu, and 
the corresponding quasiparticle spectra we find the specific heat c„(T) shown in Fig. 4 as a 
function of T. Although these parameters were chosen to fit the experimental data which is 
also shown, the agreement between theory and experiment is truly remarkable. The point to 
appreciate is that in general two different Uu and U± imply two separate transitions at Tj and 
T c x . Thus, to agree with the experiments which features a single transition at T c = 1.5 K we 
had to use both degrees of freedoms. Namely, fitting T c determined both coupling constants 
U± and U || . Now, one might suggest that it is a short-coming of our model that we have to 
rely on such an accidental coincidence of Uu and U± to fit one number T™ p . However, in the 
light of the very good fit to the very non trivial experimental variation of c v with temperature 
this is not a strong objection. Indeed, the fact that having fitted to T,f xp only and we have 
reproduced the slope as T goes to zero and the size of the jump at T c has to be regarded as 
a confirmation of our model. 



Concluding remarks. Our calculations [[15| of the superfluid density n s (T) |8| and thermal 
conductivity k(T) for the same parameters as above, provide further evidence in support 
of our model. Note that these calculations do not involve additional adjustable parameters 
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Figure 5: The structure of the order parameter. The big arrow corresponds to the orbital 
momentum I of the Cooper pair, while small arrows show the spins of the electrons which 
add to s = 1. In the frame in which l z = 1 one has s z = 0. 

and thus the results are not only qualitative consequences of nodes of the gap on the a and 
/3 sheets of the Fermi Surface but are also quantitative predictions of the theory. 

In what follows we conclude this brief survey of our very encouraging results by two further 
important comments. The first concerns the constraints A^ m , (T) = imposed at each step 
of the selfconsistency cycle during the above calculations. Without this constraint there 
would be a second transition T/ < T c , where m ,(T) becomes non zero and this would 
imply a peak in c v (T) at T[ ~ 0.2 K, which has not been seen experimentally. Fortunately, 
we have been able to show by explicit quantitative calculations that a small amount of 
disorder will eliminate the A^ m , component of the order parameters without changing the 
other amplitudes in Eq. (4) very much. The details of these calculations will be published 
elsewhere jl5j . This justifies, at the present level of the mean-free path in samples on which the 
experiments were made, the simultaneous neglect of disorder and the /-component A^ , (T). 

The other comment is that we have studied the stability of the above model (Eqs. 1-5, 
Fig. 5) to introducing further interaction constants and found it to be satisfactorily robust. 
In particular we have carried a number of calculations with ii ^ and v! 7^ and, as we 
shall report in a separate publication, we found |l5| that the overall picture presented above 
remained the same. Namely, our conclusion that a model with just two interaction constants, 
one in plane U«, and one out of plane U± (Fig. 4) roughly of the same size, around 50 meV, 
is capable of explaining all the available data we have analysed remained valid. 
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